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We simulate the bond and site percolation models on the simple-cubic lattice with linear sizes 
up to L = 512, and estimate the percolation thresholds to be pc(bond) — 0.248 81182(10) and 
Pc(site) — 0.311 607 7(2). By performing extensive simulations at these estimated critical points, we 
then estimate the critical exponents l/v = 1.1410(15), P/v — 0.477 05(15), the leading correction 
exponent yi = —1.2(2), and the shortest-path exponent dmin ~ 1.375 6(3). Various universal am- 
plitudes are also obtained, including wrapping probabilities, ratios associated with the cluster-size 
distribution, and the excess cluster number. We observe that the leading finite-size corrections in 
certain wrapping probabilities are governed by an exponent ~ —2, rather than yi « —1.2. 
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I. INTRODUCTION 

Percolation [l| is a cornerstone of the theory of critical 
phenomena 0] , and a central topic in probability 0, 0] . 
In two dimensions. Coulomb gas arguments Q and con- 
formal field theory Q predict the exact values of the bulk 
critical exponents /3 = 5/36, v = 4/3, which have been 
confirmed rigorously in the specific case of triangular- 
lattice site percolation Q. Exact values of the percola- 
tion thresholds Pc on several two-dimensional lattices are 
also known Q. In particular, it is known rigorously Q 
that Pc = 1/2 for bond percolation on the square lat- 
tice. For all d greater than or equal to the upper critical 
dimension [l^l of dc = 6, the mean- field values for the 
exponents (5 = 1, di' = 5 are believed to hold; this has 
been proved rigorously [Ullill for d > 19. 

For dimensions 2 < d < 6 by contrast, no exact values 
for either the critical exponents or percolation thresholds 
are known. Significant effort has therefore been expended 
on obtaining ever more accurate estimates, especially in 
three dimensions. 

In addition to percolation thresholds and critical ex- 
ponents, crossing probabilities [l^ [l3| also play an im- 
portant role in studies of percolation. For lattices drawn 
on a toriis, the analogous quantities are wrapping prob- 
abilities [l5|, and in two dimensions their values can be 
determined exactly [3] ■ The three-dimensional case [l^ 
has been far less studied however. Precisely estimating 
wrapping probabilities on the simple-cubic lattice repre- 
sents one of the central undertakings of the current work. 

In addition to their intrinsic importance, wrapping 
probabilities have proved to be an effective practical 
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means of estimating percolation thresholds [l^, [l^ . Us- 
ing Monte Carlo (MC) simulations and performing a 
careful finite-size scaling analysis of various wrapping 
probabilities in the neighbourhood of the transition, we 
obtain very accurate estimates of pc for both site and 
bond percolation. We observe numerically that the lead- 
ing finite-size corrections for certain wrapping probabili- 
ties appear to be governed by an exponent « —2, rather 
than by the leading irrelevant exponent yi w —1.2. 

We then estimate the thermal exponent yt = l/i^ by 
fixing p to our best estimate oipc, and studying the diver- 
gence with linear size L of the derivative of the wrapping 
probability, which is proportional to the covariance of its 
indicator with the number of bonds. We find this pro- 
cedure for estimating yt preferable to methods in which 
yt is estimated by studying how quantities behave in a 
neighbourhood of p values around Pc- In particular, we 
believe the current method produces more reliable error 
estimates. 

The remainder of this paper is organized as follows. 
The simulation method and the sampled quantities are 
discussed in Section|lTl The results for the wrapping prob- 
abilities and thresholds are given in Section Hill Critical 
exponents and the excess cluster number are discussed in 
Section Hvl We then finally conclude with a discussion in 
Section IVl 



II. SAMPLED QUANTITIES 

We study bond and site percolation on the periodic 
L X L X L simple-cubic lattice with linear system sizes 
L = 8, 12, 16, 24, 32, 48, 64, 128, 256, 512. For each 
system size, we produced at least 2.5 x 10^ independent 
samples. Each independent bond (site) configuration is 
generated by independently occupying each bond (site) 
with probability p. The clusters in each configuration are 
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TABLE I: Fits of the wrapping probabilities R^-^\ R'~°'\ R^^'' , and the ratios Qi,Q2 for bond percolation. We did not obtain 
stable fits with j/i free for i?''^'. 





-^min 




Pc 


yt 




9i 


6i 




&2 


Qi 


16 
24 
32 


53/40 
33/33 
28/26 


0.248 812 03(5) 
0.248 81198(6) 
0.248 81193(7) 


1.16(1) 
1.16(2) 
1.19(3) 


0.865 37(1) 
0.865 35(2) 
0.865 33(2) 


-0.36(1) 
-0.31(2) 
-0.31(3) 


-0.0423(5) 

-0.040(2) 

-0.036(3) 


-1.2 
-1.2 
-1.2 


0.341(5) 

0.31(2) 

0.25(5) 


16 
24 
32 


44/39 
31/32 
28/25 


0.248 81184(8) 
0.248 81188(9) 
0.248 81196(14) 


1.16(1) 
1.19(2) 
1.19(3) 


0.865 39(3) 
0.865 29(4) 
0.865 4(2) 


-0.36(1) 
-0.32(3) 
-0.31(3) 


-0.10(4) 
-0.10(8) 
-0.02(4) 


-1.34(9) 

-1.3(2) 

-1.0(5) 


0.50(8) 

0.5(2) 

0.2(3) 


Q2 


32 

48 
64 


28/25 
16/18 
10/11 


0.248 81120(5) 
0.248 81195(6) 
0.248 81184(11) 


1.17(2) 
1.14(2) 
1.12(3) 


0.633 58(3) 
0.633 50(3) 
0.633 4(2) 


-0.80(5) 
-0.89(8) 
-1.0(2) 


-0.104(4) 
-0.088(9) 
-0.05(4) 


-1.2 
-1.2 
-1.2 


0.05(7) 
-0.3(2) 
-1(1) 


32 
48 
64 


28/26 
16/19 
10/12 


0.248 812 02(6) 
0.248 81193(7) 
0.248 81182(11) 


1.17(2) 
1.14(2) 
1.12(3) 


0.633 58(5) 
0.633 46(7) 
0.633 3(2) 


-0.80(5) 
-0.89(8) 
-1.0(2) 


-0.097(8) 

-0.15(4) 

-0.5(6) 


-1.08(3) 
-1.22(7) 
-1.5(4) 






16 

24 
32 


41/37 
30/31 
25/24 


0.248 81181(4) 
0.248 81183(4) 
0.248 81182(6) 


1.143(7) 

1.15(2) 

1.15(2) 


0.257 77(2) 
0.257 78(3) 
0.257 76(5) 


-1.22(3) 
-1.22(6) 
-1.20(8) 


0.005(2) 
0.003(3) 
0.006(7) 


-1.2 
-1.2 
-1.2 


-0.23(1) 
-0.26(4) 
-0.20(10) 


16 
24 
32 


41/37 
31/31 
25/24 


0.248 81182(4) 
0.248 81184(4) 
0.248 81182(6) 


1.144(7) 

1.15(2) 

1.15(2) 


0.257 79(2) 
0.257 79(2) 
0.257 77(4) 


-1.22(3) 
-1.22(6) 
-1.20(8) 


0.18(2) 
0.22(8) 
0.1(1) 


-1.83(4) 

-1.9(2) 

-1.7(3) 






16 

24 
32 


40/39 
25/32 
22/25 


0.248 81182(4) 
0.248 81182(5) 
0.248 81183(6) 


1.149(7) 

1.14(2) 

1.14(2) 


0.459 99(3) 
0.459 97(5) 
0.459 98(7) 


-1.65(4) 
-1.74(9) 
-1.7(2) 


0.004(2) 
0.003(4) 
0.005(9) 


-1.2 
-1.2 
-1.2 


0.73(2) 
0.72(6) 
0.7(2) 


16 
24 
32 


40/39 
25/32 
22/25 


0.248 81182(4) 
0.248 81182(4) 
0.248 81182(5) 


1.149(7) 

1.14(2) 

1.14(2) 


0.459 97(2) 
0.459 95(3) 
0.459 96(5) 


-1.65(4) 
-1.74(9) 
-1.74(2) 


0.81(6) 

0.8(2) 

1.0(9) 


-2.06(3) 
-2.05(8) 
-2.1(3) 






16 

24 
32 


44/38 
35/31 
23/24 


0.248 81185(6 ) 
0.248 81191(6 ) 
0.248 81185(8 ) 


1.14(1) 
1.15(2) 
1.17(3) 


0.080 41(2) 
0.080 43(3) 
0.080 39(4) 


-0.66(2) 
-0.63(5) 
-0.59(5) 


0.010(1) 
0.007(3) 
0.014(6) 


-1.2 
-1.2 
-1.2 


-0.076(8) 

-0.04(3) 

-0.15(9) 



identified using breadth-first search. The number of sites 
in each cluster defines its size. 

We measured the following observables in our simula- 
tions: 

• The number of occupied bonds J\fb for bond perco- 
lation, and the number of occupied sites Afs for site 
percolation; 

• The number of clusters A/'c; 

• The size Ci of the largest cluster; 

• The cluster-size moments Sm ~ J2c IC"!™ with 
m = 0, 2, 4. The sum runs over all clusters C, and 
So is simply the number of clusters; 

• An observable S := max max d(xc,v) used to de- 

c yec ^ 

termine the shortest-path exponent. Here d{x, y) 
denotes the graph distance from site x to site y, and 
xc is the vertex in cluster C with the smallest ver- 
tex label, according to some fixed (but arbitrary) 
vertex labelling. 

• The indicators 7^(^^ 7^(^'^ and 7^(^), for the event 
that a cluster wraps around the lattice in the x, y, 
or z directions, respectively. 



From these observables we calculated the following 
quantities: 

• The mean size of the largest cluster Ci = (Ci), 
which at pc scales like Ci ^ L^'' with yh = df = 
d — where d/ is the fractal dimension; 

• The cluster density p = {Mc)/L'^] 

• The mean size of the cluster at the origin, x = 
{S2)/L''-, which at pc scales like x ~ L'^Vh-d.-^ 

• The dimensionless ratios 



'^i ^ /n \9 ' — 777771 



(352^ - 254) 



(1) 



• The shortest-path length S = (5), which at pc 
scales like S ~ /^t^min -^^ith (i,„in the shortest-path 
fractal dimension; 

• The wrapping probabilities 

R^"'^ =1 - ((1 - 7e("')(i - n'^y^){i - n'-^^)) , (2) 
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Here gives the probability that a winding 

exists in the x direction, i?^"^ gives the proba- 
bihty that a winding exists in at least one of 
the three possible directions, and R^^^ gives the 
probability that windings simultaneously exist in 
all three possible directions. Near pc, we expect 
each of these wrapping probabilities to behave as 
~ f{{p — pc)Ly*), where / is a scaling function; 



The covariance of 7^^^' and A/j, 

g(^) = (7^(-)A/■b)-(7^(-))(AA,) 



= p{l -p) 



dp 



(3) 



(x) 

At Pc, we expect ^ L^*. An analogous defini- 
tion of g^^' , with J\f{, being replaced with Afg , was 
used for site percolation. 

To derive ([3]), one can explicitly differentiate (T?,'^^') 
with respect to p, and use the fact that (A/t,) = p\E\ 
where \E\ is the total number of edges on the lattice. 

The complete set of data for all obscrvablcs, for 
both bond and site percolation, is contained as a file 
percolation . ta r . g z in the preprint version of this pa- 
per at |arXiv.org| |28| . 



III. ESTIMATING pc 
A. Bond percolation 

We estimate the thresholds of bond and site percola- 
tion by studying the finite-size scaling of the wrapping 
probabilities R^'^\ R^°'\ R^^\ and the dimensionless ra- 
tios Qi, Q2- Around Pc, we perform least-squares fits of 
the MC data for these quantities by the ansatz 



(4) 



where e = pc—p, Oc 'is a universal constant, and yi is the 
leading correction exponent. We perform fits with both 
bi and 62 free, as well as fits with 62 being set identically 
to zero. By performing fits with yi free we estimate that 
yi = —1.2(2). We also perform fits with yi fixed to y,; = 
-1.2. 

As a precaution against correction-to-scaling terms 
that we have neglected in our chosen ansatz, we impose 
a lower cuttoff L > L^un on the data points admitted 
in the fit, and we systematically study the effect on the 
value of increasing Lmin- In general, our preferred fit 
for any given ansatz corresponds to the smallest L^in for 
which divided by the number of degrees of freedom 
(DF) is 0(1), and for which subsequent increases in imin 
do not cause to drop by vastly more than one unit per 
degree of freedom. 



0.259 



^0.258 
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FIG. 1: Plots of 7?''^'(p,L) (top) and R'-''\p,L) (bottom) vs 
L for fixed values of p, for bond percolation. In both cases, 
the curves correspond to our preferred fit of the MC data for 
R{p, L) by the ansatz @ ; the dashed curve corresponds to set- 
ting p = 0.248 811 82. The blue strips indicate an interval of 
one sigma above and below the estimates 7?i^^ = 0.257 78(6) 



and R, 



(a) 



0.459 97(8). 



Table U summarizes the results of these fits. From the 
fits, we can see the finite-size corrections of Qi and Q2 
are dominated by the exponent yi ~ —1.2. From Qi and 
Q2, we estimate Pc = 0.248 811 9(3), and their universal 
critical values = 0.865 4(2) and Q2,c = 0.633 5(2). 

For and R^'^\ fixing yi = — 1.2 and including both 
the bi and 62 terms we find that 61 is consistent with 
zero, while 62 is clearly nonzero. Furthermore, if we set 
62 = and leave yi free, we find yi « —2. This suggests 
that either the amplitudes of the leading corrections of 
and vanish identically, or at least that they 
are sufficiently small that they cannot be detected from 
our data. Due to these weak finite-size corrections, the 
values of Pc fitted from i?^^) and i?*^"-* are much more 
stable than those obtained from Qi and Q2- From 
and we estimate Pc = 0.248 81182(10). For R'^^\ 

we report only the fits with corrections biL~^-^ + 62^"^. 
If yi is left free the fits become unstable, regardless of 
whether the b2L~^ term is included. From R^'^\ we esti- 
mate Pc = 0.248 81185(15) which is consistent with the 
value obtained from i?*^^^ and i?*^") . From these fits, we es- 



timate the universal wrapping probabilities to be i?, 



(x) 



0.257 78(6), = 0.459 97(8) and R'^c'' = 0.080 41(8) 



(3) 
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TABLE II: Fits of the wrapping probabilities R^^\ R^°'\ R^''^'' , and the ratios Q\,Q2 for site percolation. For _R'^' we obtain 
unstable results when yi is free. 





-^min 


XVDF 


Pc 


yt 




11 


bi 


Vi 


h2 




32 


19/16 


0.311606 9(2) 


1.14(2) 


0.865 05(2) 


-0.22(2) 


0.062(2) 


-1.2 


0(3) 




48 


11/11 


0.3116070(2) 


1.11(3) 


0.865 09(3) 


-0.25(3) 


0.054(6) 


-1.2 


0.2(2) 




64 


3/6 


0.311607 7(3) 


1.12(6) 


0.865 26(7) 


-0.24(6) 


0.01(2) 


-1.2 


1.4(5) 




32 


19/16 


0.311606 9(2) 


1.15(2) 


0.865 06(3) 


-0.22(2) 


0.063(4) 


-1.11(2) 


- 




48 


10/11 


0.3116071(2) 


1.11(3) 


0.865 12(4) 


-0.25(3) 


0.09(2) 


-1.22(7) 


- 




64 


3/6 


0.311607 7(3) 


1.12(6) 


0.865 27(5) 


-0.24(6) 


0.9(10) 


-1.8(3) 


- 




64 


3/6 


0.3116076(2) 


1.12(4) 


0.633 3(1) 


-0.56(9) 


0.02(3) 


-1.2 


5.1(7) 


02 


48 


13/11 


0.3116072(1) 


1.14(2) 


0.633 06(4) 


-0.52(4) 


0.9(1) 


-1.52(3) 


- 




64 


2/6 


0.3116076(2) 


1.12(4) 


0.633 29(8) 


-0.56(9) 


4(2) 


-1.9(2) 


- 




16 


42/39 


0.311607 85(5) 


1.13(1) 


0.257 89(2) 


-0.76(4) 


0.004(1) 


-1.2 


-0.22(1) 


^(-) 


24 


30/31 


0.311607 74(6) 


1.14(2) 


0.257 84(3) 


-0.75(5) 


0.009(2) 


-1.2 


-0.29(3) 




32 


24/24 


0.311607 66(7) 


1.14(2) 


0.257 80(3) 


-0.73(5) 


0.015(4) 


-1.2 


-0.39(6) 




16 


39/40 


0.311607 70(5) 


1.12(2) 


0.460 02(2) 


-1.09(6) 


0.023(2) 


-1.2 


0.08(2) 




24 


25/32 


0.311607 67(7) 


1.13(2) 


0.459 99(4) 


-1.05(6) 


0.025(3) 


-1.2 


0.05(4) 




32 


19/24 


0.311607 65(8) 


1.13(2) 


0.459 98(5) 


-1.06(7) 


0.027(6) 


-1.2 


0.02(9) 




16 


36/40 


0.311607 75(6) 


1.12(2) 


0.460 06(3) 


-1.09(6) 


0.055(5) 


-1.33(4) 






24 


25/32 


0.311607 68(8) 


1.13(2) 


0.460 01(5) 


-1.05(7) 


0.039(9) 


-1.21(9) 






32 


19/24 


0.311607 65(9) 


1.13(2) 


0.459 99(7) 


-1.06(7) 


0.03(2) 


-1.1(2) 






16 


50/38 


0.311608 01(8) 


1.14(2) 


0.080 55(1) 


-0.38(3) 


-0.010(8) 


-1.2 


-0.30(1) 




24 


27/30 


0.311607 79(9) 


1.14(2) 


0.080 49(2) 


-0.39(4) 


-0.004(2) 


-1.2 


-0.38(3) 




32 


18/23 


0.311607 65(11) 


1.15(3) 


0.080 45(3) 


-0.38(4) 


-0.002(3) 


-1.2 


-0.47(5) 


16 


40/38 


0.311607 89(7) 


1.15(2) 


0.080 510(9) 


-0.38(3) 


-0.21(1) 


-1.77(2) 






24 


26/30 


0.311607 77(8) 


1.14(2) 


0.080 48(2) 


-0.39(4) 


-0.30(5) 


-1.88(5) 






32 


18/23 


0.311607 66(10) 


1.15(3) 


0.080 46(2) 


-0.38(4) 


-0.6(2) 


-2.1(2) 





In Fig. [1] we illustrate our estimate of pc by plotting 
i?*^^) and i?^"^ vs L. Precisely ai p = p^, as, L 00 
the data should tend to a horizontal line, whereas the 
data with p ^ Pc will bend upward or downward. Fig. [1] 
suggests the central value 0.248 812 6 from [23| is an over- 
estimate. 

In Fig. H we plot the data at p = 0.248 811 8 for 
and i?'"-' vs and for Qi and Q2 vs L~^ '^ . The figure 
strongly suggests that the correction L~^-^ dominates in 
Qi and (52, but vanishes (or is very weak) in R'^^^ and 



B. Site percolation 

For site percolation, we again estimate pc by fitting 
Qi, Q2, and by Eq. Q. The fitting 

procedure is similar to that of bond percolation, and the 
results are summarized in Table |lll From the table, we 
can see the fits of Qi and Q2 are less stable for site perco- 
lation than bond percolation. The ratio x^/DF remains 
large until imin > 32 for Qi and Lmin > 48 for Q2, and 
the resulting estimates of Pc range from 0.311 606 9(2) to 
0.311607 7(3). 

The fits of the wrapping probabilities are better be- 
haved, as was the case for bond percolation. For R^'^\ 
fixing yi — —1.2 and including both the 61 and 62 terms, 
we find that hi is consistent with zero, while 62 is clearly 
nonzero. Furthermore, if we set 62 = and leave y,; free. 



we find Hi ~ —2. This suggests that the amplitude of the 
leading correction of R^'^^ is smaller than the resolution 
of our fits, and might possibly be zero. The fits of the 
/J^"^ data, however, quite clearly indicate the presence of 
the biL~^-^ term. For we report only the fits with 

corrections biL~^-^ + b2L~^; if yi is left free the fits be- 
come unstable, regardless of whether the 62^"^ term is 
included. As for R^°'\ the amplitude 61 appears to take 
a nonzero value. These observations suggest the leading 
correction L~^'^ does not gencrically vanish for all wrap- 
ping probabilities, but rather that the amplitudes in some 
cases are smaller than the resolution of our simulations. 

Comparing the various fits, we estimate Pc = 
0.311607 7(2) for site percolation, which is consistent 
with the previous result 0.311607 7(4) [Hj. In ad- 
dition, we estimate the universal wrapping probabili- 
ties to be Ri""^ = 0.25782(6), = 0.459 99(8), and 
= 0.080 46(6), which are consistent with those esti- 
mated from bond percolation. In Fig. [31 we show plots 
of i?^^-* and i?'-"'' which illustrate our estimate of pc- 



IV. RESULTS AT pc 

In this section, we estimate the critical exponents yt, 
yh, and dmin, as well as the excess cluster number. Fixing 
p at our estimated thresholds for bond and site percola- 
tion, we study the covariances and (?[^\ the mean 
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0.002 0.004 0.006 

FIG. 2: Plots of Qi and Qa vs L'^'^ (top), and and Ti^"' 
vs (bottom), with p — 0.248 811 8, for bond percolation. 
The solid lines are simply to guide the eye. 

size of the largest cluster Ci , the mean size of the clus- 
ter at the origin x, the shortest-path length S, and the 

(a:) fee) 

cluster density p. The MC data of g^^, , Ci, x and 
S are fitted by the ansatz 

A = Ly^{ao + hiL-^-'^ +b2L-'^) . (5) 

We perform fits using different combinations of the two 
corrections 6iL~^'^ and 62^"^ and compare the results. 

A. Estimating yt 

(x) 

We estimate yt by studying the covariances and 

9sR^ both of which scale as ~ L^* at the critical point. 
We find this procedure for estimating yt preferable to 
methods, such as that employed in [2lj, in which yt is es- 
timated by studying how quantities behave in the neigh- 
bourhood of Pc as the system deviates from criticality. In 
particular, we believe the current method produces more 
reliable error estimates. 

We fit the data for 5^^^ at p = 0.248 8118 and 5^^^ 
at p = 0.311 607 7 to Eq. (O, and the results are shown 
in Table IIIII The estimate of yt from g^^j^ produces a 
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FIG. 3: Plots of i?'^)(p,L) (top) and R^^^p, L) (bottom) 
vs L for fixed values of p, for site percolation. In both cases, 
the curves correspond to our preferred fit of the MC data for 
R{p, L) by ansatz Q; the dashed curve corresponds to setting 
p = 0.311607 7. The blue strips indicate an interval of one 
sigma above and below the estimates R'f ^ = 0.257 82(6) and 
Ri"'^ = 0.459 99(8). 



TABLE III: Fits of covariances and 









yt 


ao 


hi 


fc2 




16 


4/4 


1.140 4(9) 


0.231(1) 


-0.03(2) 


0.1(2) 


(^) 


24 


4/3 


1.140 6(13) 


0.231(2) 


-0.02(5) 


0.0(4) 


9lii 


16 


4/5 


1.140 9(4) 


0.230 7(3) 


-0.012(3) 






24 


4/4 


1.140 6(6) 


0.231 1(6) 


-0.017(7) 






16 


5/4 


1.1416(4) 


0.1551(3) 


-0.004(7) 


-0.06(5) 


(i) 


24 


4/3 


1.141 1(6) 


0.155 4(6) 


-0.02(2) 


-0.1(2) 




16 


7/5 


1.141 1(2) 


0.155 5(1) 


-0.013(1) 






24 


4/4 


1.1414(3) 


0.155 3(2) 


-0.010(2) 





smaller error bar than that from . From these fits 
we take our final, somewhat conservative, estimate to be 
yt = 1.1410(15). 

In Fig. m we plot (In ,g[^^ - yt In L) and (In g'^^j^ - yt In L) 
vs InL using three different values of yt'. our estimate, 
as well as our estimate plus or minus three standard de- 
viations. Using the true value of yt should produce a 
horizontal line for large L. In the figure, the data using 
yt — 1.136 5 and yt = 1.145 5 respectively bend upward 
and downward, suggesting that the true value of yt does 
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FIG. 4: Plots of (In p^^^ - yt In L) (top) and (In - yt In L) 
(bottom) vs InL illustrating our estimate yt = 1.1410(15). 
The dashed curves are simply to guide the eye. 



FIG. 5; Plots of (InCj - yh InL) (top) and (InCf - yu InL) 
(bottom) vs InL to show our estimate yh ~ 2.522 95(15). The 
dashed curves are simply to guide the eye. 



TABLE IV: Fits of Ci and x- The superscripts b and s denote 
bond and site percolation, respectively. 







XVDF 


Vh 


0,0 


bi 


62 




16 
24 


3/4 
3/3 


2.522 86(5) 
2.522 89(7) 


0.9394(3) 
0.939 3(4) 


-0.014(6) 
-0.009(11) 


0.22(4) 
0.2(1) 


24 
32 


5/4 
3/3 


2.522 98(3) 
2.522 94(4) 


0.938 8(2) 
0.939 0(2) 


0.009(2) 
0.005(3) 






16 
24 


4/4 
3/3 


2.523 03(4) 
2.523 00(5) 


1.125 7(5) 

1.126 2(7) 


0.14(1) 
0.12(2) 


0.18(7) 
0.3(2) 


24 
32 


6/4 
4/3 


2.523 08(3) 
2.523 05(3) 


1.1251(3) 
1.125 5(4) 


0.157(4) 
0.151(6) 




CI 


16 
24 


5/4 
5/3 


2.522 99(3) 

2.523 00(5) 


0.471 16(7) 
0.471 1(2) 


0.024(2) 
0.024(4) 


-0.44(2) 
-0.45(4) 




32 
48 


0.9/2 
0.7/1 


2.522 91(5) 
2.522 94(9) 


0.2841(2) 
0.284 0(4) 


-0.001(7) 
-0.007(18) 


-1.15(9) 
-1.3(3) 


32 
48 


0.9/3 
0.9/2 


2.522 92(1) 
2.522 91(2) 


0.284 06(3) 
0.284 08(7) 




-1.16(1) 
-1.17(5) 



indeed lie within 3 sigmas of our estimate. The data 
with yt = 1.141 appear to be consistent with an asymp- 
totically horizontal line. We note that while the curve 
appears to be increasing around the point at L = 512 for 
bond percolation, it instead slightly decreases for site per- 
colation, suggesting that in fact this movement is domi- 
nated (or even entirely caused) by noise. 



B. Estimating yh 

We estimate yh by studying the divergence of Ci and 
X as L increases with p fixed to our best estimates of 
Pc. Wc fit the MC data for Ci and x by Eq. (0, with 
the exponent then corresponding to yh and 2yh — 
d, respectively. The results are reported in Table IIVI 
We use superscripts b and s to distinguish bond and site 
percolation. For and x^i the amplitude 61 is quite 
small, while 61 in and Cf is clearly present. In the fits 
of X* with one correction term biL~^-^, the ratio x^/DF 
remains large until Lmin > 64. We therefore show the 
fits with correction 62^"^ instead. Comparing these fits, 
wc estimate yh = 2.522 95(15). 

In Fig. [51 we plot (In — yh In L) and (In Cf — yh In L) 
vs InL using three different values of yh- our estimate, 
as well as our estimate plus or minus three standard de- 
viations. As L increases, the data with yh = 2.522 50 
and 2.523 40 respectively slope upward and downward, 
while the data with yh = 2.522 95 are consistent with an 
asymptotically horizontal line. 



C. Estimating dmin 

We estimate the shortest-path fractal dimension dmin 
by studying the quantity S at our estimated thresholds. 
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TABLE V: Fits of S. The superscripts b and s denote bond 
and site percolation, respectively. 







X7DF 




ao 


bi 


62 


s" 


24 


2/3 


1.375 26(5) 


1.814 9(5) 


-0.65(2) 


-3.8(2) 




32 


1/2 


1.375 33(7) 


1.814 2(7) 


-0.59(5) 


-4.4(4) 




48 


0/2 


1.375 30(9) 


1.815(1) 


-0.63(9) 


-4(1) 




16 


5/4 


1.375 80(2) 


1.383 4(2) 


-3.432(5) 


2.72(3) 




24 


4/4 


1.375 77(3) 


1.383 6(3) 


-3.45(2) 


2.82(3) 




32 


4/2 


1.375 76(5) 


1.383 7(4) 


-3.45(3) 


2.9(3) 



TABLE VL Fits of p. The superscripts b and s denote bond 
and site percolation, respectively. 







X7DF 


Pc 


6 


bi 




16 


3/5 


0.272 932 83(1) 


0.679(3) 


0.1(6) 


p' 


24 


1/4 


0.272 932 83(1) 


0.674(6) 


3(4) 




16 


2/7 


0.272 932 83(1) 


0.678 9(6) 






24 


2/6 


0.272 932 83(1) 


0.679(2) 






12 


4/6 


0.052 438 218(3) 


0.674 5(5) 


0.02(8) 




16 


4/5 


0.052 438 218(3) 


0.674 7(8) 


-0.02(21) 




24 


4/4 


0.052 438 218(3) 


0.674(2) 


0.2(10) 




12 


4/7 


0.052 438 218(3) 


0.674 6(2) 






16 


4/6 


0.052 438 218(3) 


0.674 6(3) 






24 


4/5 


0.052 438 218(3) 


0.674 6(5) 





The MC data for S are fitted to Eq. ^ with the expo- 
nent yj{ replaced by dmin, and the results are reported 
in Table |Vl We again use the superscripts b and s to 
distinguish bond and site percolation. In the fits, both 
bi and 62 are clearly observable for and S*^. And 
when we set 62 = 0, the ratio /DF remains relatively 
large. We also did the fits by replacing the correction 
with &2 by a constant term cq in Eq. ([5]), and obtained 
dmin(bond) = 1.375 55(6) and dnui„(site) = 1.375 59(6). 
Comparing these fits, we estimate dmin = 1.375 6(3). 

To illustrate this estimate, Fig. |6] shows a log-log plot 
of S versus L. 
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FIG. 6: Log-log plot of 5* versus L for bond and site perco- 
lation. Two straight lines with slope 1.375 6 are included for 
comparison. 



D. Excess number of clusters 



The cluster density tends to a finite limit pc = 
limi_j.oo limp-!-pc P at criticality. While the value of 
Pc is non-universal, the excess cluster number b := 
lim^^oo limp-^pc L'^{p~ Pc) is universal [l^]. To estimate 
b, we study p with p fixed to our estimated thresholds for 
bond and site percolation and fit the data to the ansatz 



L- 



biL-^) 



(6) 



The resulting fits are summarized in Table I VII where 
we again use superscripts b and s to differentiate the 
bond and site cases. We report fits both with bi free 
and &i = 0. We find that p can be well fitted to ^ 
with 61 = fixed. Leaving 61 free, we find that bi is 
consistent with zero, suggesting that the leading correc- 
tion exponent might be even smaller than —2. We also 
performed fits in which the leading correction exponent 
was fixed to —1.2 and —3, and in both cases the resulting 
estimates of pc and b were consistent with those reported 
in Table rvTl Leaving the leading correction exponent free 
produces unstable fits however. Comparing these fits, we 
estimate b = 0.675(2). 

To our knowledge, no estimates of b on the periodic 
L X L X L simple cubic lattice have previously been re- 
ported in the literature; on the L x L square lattice 
b = 0.883 5(8) [13 ■ The excess cluster number was stud- 
ied in on an L X L X L' lattice with L' ^ L. Naively, 
extrapolating their results to L' = L gives an estimate 
of 6 « 0.412 which is significantly below our estimate. 
We also note that our estimate of the number of clus- 
ters Pc = 0.272 932 83(1) differs slightly with the estimate 
Pc = 0.272 9310(5) reported in 



V. DISCUSSION 

We study in this paper standard bond and site perco- 
lation on the three-dimensional simple-cubic lattice with 
periodic boundary conditions. Using extensive Monte 
Carlo simulations and finite-size scaling analysis, we re- 
port the estimates: Pc = 0.248 81182(10) (bond) and 
Pc = 0.311607 7(2) (site). The bulk thermal and mag- 
netic exponents are estimated to be yt = 1.141 0(15) and 
yh = 2.522 95(15), the shortest-path fractal dimension 
to be dniin = 1.375 6(3), and the leading irrelevant ex- 
ponent to be yi = —1.2(2). The universal value of the 
excess cluster number is estimated to be 6 = 0.675(2). 

We emphasize that the reported estimates of pc are 
obtained by studying wrapping probabilities, which are 
found to have weaker corrections to scaling than dimen- 
sionless ratios constructed from moments of magnetic 
quantities such as Ci and Sm- In particular, we find evi- 
dence suggesting the leading correction exponent in cer- 
tain wrapping probabilities (i?*^^^ and for bond per- 
colation, for site percolation) may be w — 2 rather 
than —1.2, although the reasons are not clear. The uni- 
versal values of the wrapping probabilities we studied are 
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TABLE VII: Summary of estimated thresholds, critical exponents, universal wrapping probabilities, and excess cluster number 
of bond and site percolation on the simple-cubic lattice. 

Rof. pclbond) pTCsito) = \ j v yu = dj d~ 'W^ h 

7201 0.2488126(5) 1.12(2) 2.523(4) 

[221 0.3116080(4) 

p3l 0.3116081(13) 1.141(2) 2.5230(3) -1.61(13) 

[171 0.249 0(2) 0.3115(3) 1.15(2) 0.265(6) 0.471(8) 0.084(4) 

pTl 0.3116077(4) 1.1450(7)° 2.5226(1)° 

[241 0.248 812 0(5) 1.142(3) 2.523 5(8) 

[251 1.375 6(6) 

[261 1.142(8) -1-0(2) 

This work 0.248 81182(10) 0.3116077(2) 1.1410(15) 2.52295(15) 1.3756(3) -1.2(2) 0.25780(6) 0.45998(8) 0.08044(8) 0.675(2) 

"The values of yt and j/^ reported in the abstract of [2l]| contained 
typographical errors. 



estimated to be: Rc^ ^ 0.257 80(6), RT' 0.459 98(8), 
and Rc = 0.080 44(8), by comparing the results for 
bond and site percolation. 

From these values we can estimate other wrapping 
probabilities discussed in the literature, such as 

: = (7^(^)(l_7^(^))(l-7^("))) 

= i(2i?(") -3i?(^)) , 

3 

= ^(3i?(^) -2i?(3) ^ 

: = (7^(^'7^(■^)) = i(3i?(^' + - i?(°^) . 

o 

In words, i?*-^^ is the probability that a winding exists in 
one given direction but not in the other two directions; 
i?*-^' is the probability that a winding exists in two given 
directions but not in the third; and R^^'V^ is the probabil- 
ity that a winding exists in two given directions, regard- 
less of whether a winding exists in the third direction. 



We obtain R^c^ = 0.075 67(14), R^c^ = 0.050 85(14), and 
R^^^y^ = 0.13129(12). 

Table IVIII summarizes the estimates presented in this 
work. For comparison, we also provide an (incomplete) 
summary of previous estimates. 
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